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Abstract Magnetic fields likely play a key role in the 
dynamics and evolution of protoplanetary disks. They 
have the potential to efficiently transport angular mo- 
nicntuin by MHD turbulence or via the magnetocen- 
trifugal acceleration of outflows from the disk surface. 
Magnetically-driven mixing has implications for disk 
chemistry and evolution of the grain population, and 
the effective viscous response of the disk determines 
whether planets migrate inwards or outwards. How- 
ever, the weak ionisation of protoplanetary disks means 
that magnetic fields may not be able to effectively cou- 
ple to the matter. I examine the magnetic diffusivity 
in a minimum solar nebula model and present calcula- 
tions of the ionisation equilibrium and magnetic diffu- 
sivity as a function of height from the disk midplane at 
radii of 1 and 5 AU. Dust grains tend to suppress mag- 
netic coupling by soaking up electrons and ions from 
the gas phase and reducing the conductivity of the gas 
by many orders of magnitude. However, once grains 
have grown to a few microns in size their effect starts 
to wane and magnetic fields can begin to couple to the 
gas even at the disk midplane. Because ions are gener- 
ally decoupled from the magnetic field by neutral col- 
lisions while electrons are not, the Hall effect tends to 
dominate the diffusion of the magnetic field when it is 
able to partially couple to the gas, except at the disk 
surfaces where the low density of neutrals permits the 
ions to remain attached to the field lines. 

For a standard population of 0.1/im grains the ac- 
tive surface layers have a combined column Eactive ~ 
2gcm~^ at lAU; by the time grains have aggregated 
to 3 /im, Eactivo ~ 80gcm~^. Ionisation in the active 
layers is dominated by stellar x-rays. In the absence 
of grains, x-rays maintain magnetic coupling to 10% of 
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the disk material at lAU (i.e. Sactive ~ 150gcm~^). 
At 5AU Sactive ~ ^totai once grains have aggregated 
to 1 /Ltm in size. 

Keywords accretion, accretion disks - MHD - molec- 
ular processes - stars: formation 

1 Introduction 

Magnetic fields may efficiently transport angular mo- 
mentum in protoplanetary disks either via MHD turbu- 
lence driven by the magnetorotational instabiUty (Bal- 
bus & Hawley 1991; Hawlcy, Gammie & Balbus 1994), 
by accelerating a wind from the disk surfaces (Bland- 
ford & Payne 1982; Wardle & Konigl 1993), or by loop- 
ing above the disk surface and linking different radii. 
Advection and/or stirring by magnetically-driven tur- 
bulence affects disk chemistry (e.g. Semenov, Wiebe 
& Henning 2006; Ilgner & Nelson 2006), the aggrega- 
tion and settling of small particles that are the prelimi- 
nary stages of planet building (Johansen & Klahr 2005; 
Turner et al. 2006; Promang & Papaloizou 2006; Ciesla 
2007), and magnetic activity at the disk surface may 
produce a corona (e.g. Fleming & Stone 2003) and af- 
fect observational signatures of protoplanetary disks. A 
magnetically-mediated effective viscosity would also de- 
termine the response of the disk to the gravitational dis- 
turbances of forming planets, modifying the balance be- 
tween the torques applied by the outer and inner neigh- 
bouring regions of the disk and therefore the rate and 
direction of planetary migration (Matsumura & Pudritz 
2003; Johnson, Goodman & Menou 2006; Chambers 
2006). Finally, magnetic fields are likely to play a role 
in giant planet formation and the formation and evolu- 
tion of their satellite systems (Quillen & Trilling 1998; 
Fendt 2003). 
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One critical uncertainty is the magnitude and nature 
of the coupling between magnetic fields and protoplan- 
etary disks, which are very weakly ionised on account of 
the self-shielding against ionising sources provided by 
their high column density, and by rapid recombinations 
because of their high number densities. In addition the 
mobility of charged particles is reduced by the rela- 
tively high density, and dust grains sweep up charges 
and render them effectively immobile. Indeed, mag- 
netic fields were expected to be negligible in protoplan- 
etary disks, as extremely high diffusion was thought to 
eliminate the magnetic gradients responsible for mag- 
netic stresses except at the disk surfaces (e.g. Hayashi 
1981). But these expectations are too pessimistic. At 
any given radius in the disk the density declines rapidly 
with height because of tidal squeezing by the gravita- 
tional field of the central star. The corresponding rapid 
decline in diffusivity towards the disk surface is aided by 
external ionisation by cosmic rays and by stellar x-rays 
which dominate cosmic rays by five orders of magnitude 
within a few gcm^^ of the disk surface (Glassgold, Na- 
jita & Igea 1997; Igea & Glassgold 1999), so that at least 
the surface layers may be magnetically active (Gam- 
mie 1996; Wardle 1997). In addition, the Hall effect 
provides a dissipationless diffusion pathway that can 
maintain field gradients under a much broader range of 
conditions than would otherwise occur (Wardle & Ng 
1999; Wardle 1999; Balbus & Terquem 2001). Finally, 
aggregation and settling of grains increases the number 
and mobility of charged species in the gas phase and 
reduces diffusivity to the point that magnetic stresses 
are important even at the disk midplane (Sano & Stone 
2002). 

Here I explore the nature and magnitude of magnetic 
diffusion in protoplanetary disks through calculations of 
the ionisation equilibrium in a stratified disk exposed 
to cosmic rays and to x-rays from the central star. 



2 Minimum-mass solar nebula 

To estimate the degree of coupling of magnetic fields to 
the matter in protoplanetary disks we first need a no- 
tional model for the physical conditions within the disk. 
I adopt the minimum-mass solar nebula model (Wei- 
denschilling 1977; Hayashi 1981), in which the surface 
density of the solar nebula, S, is estimated by adding 
sufficient hydrogen and helium to the solid bodies in 
the solar system to recover standard interstellar abun- 
dances and spreading this material smoothly in a disk. 
This yields an estimate of the minimum surface density 
needed in the solar nebula to form the solar system. 



1700 r^u^' gcm-^ 



where tau is the radial distance r from the Sun in as- 
tronomical units. The disk mass enclosed within r is 



Md = 27rr'S; dr' w 0.024 



100 AU 



1/2 



Mr. 



(2) 



where r is assumed to be much greater than the inner 
disk radius rmin. The temperature profile of the disk 
is estimated by considering the thermal balance for a 
black body of radius a at a distance r from the Sun, 



na 



Ana (JsbT , 



(3) 



where ctsb is the Stefan-Boltzmann constant, which 
yields 



T w 280 r^y^ K. 



(4) 



These rough estimates of the distribution of matter 
and temperature are broadly consistent with continuum 
observations of protoplanetary disks around low mass 
stars, albeit with resolutions equivalent to ^ 100 AU 
scales (e.g. Kitamura et al. 2002; Andrews & Williams 
2005). 

The isothermal sound speed in the disk implied by 
this temperature profile is 



Cs w 0.99 r^y'* kms^i 



(5) 



which is small compared to the local Keplerian speed 
vk, so the disk is thin (see eq [S] below). Here I have 
adopted a mean molecular weight of 7/3 appropri- 
ate for a mixture of 80% H2 and 20% He by number. 
Toomre's Q parameter (Toomre 1964) is 



ttGE 



56 r 



^1/4 

AU • 



(6) 



where f2 — vk /r is the local Keplerian frequency and G 
is the gravitational constant, so the disk is not unsta- 
ble to self gravity (as expected given Md <^ M0). Of 
course this represents the minimum solar nebula; the 
surface density could be increased by factor of ^--^30 be- 
tween 1-100 AU before the disk becomes gravitationally 
unstable (e.g. Durisen et al. 2007), although as noted 
above large disk masses do not seem to be present in 
low mass YSOs. 

Assuming that the temperature does not depend on 
height z above the disk midplane, the balance between 
tidal squeezing by the central object and its internal 
pressure imply that the gas density is stratified as 



p(r, z) = pQ{r) exp 



(1) 



2/i2 



(7) 
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where the disk scale height h is given by 

h _ 
r VK 



0.03 r'/^ . 



(8) 



The density at the midplane, pQ, is determined by eq 
([l]) and the relation S — V2TrpQh, yielding a number 
density of hydrogen nuclei 



nu = 5.8 X lO^"* r^u cm"^ 



(9) 



at the midplane. It should be emphasised that the tem- 
perature certainly depends on height because of the ef- 
fects of external irradiation, height-dependent dissipa- 
tion of accretion energy and the effects of chemistry and 
optical depth on the ability of a gas parcel to cool (see 
e.g. Dullemond et al. 2007). Nevertheless the decline in 
pressure with height required by tidal confinement of 
the disk is provided primarily by the drop in density. 



3 Magnetic field strength 

There are as yet, no direct measurements of the field 
strength in protoplanetary disks, apart from remanent 
magnetism in meteorites which indicate field strengths 
in the solar nebula of 0.1-1 G at 1 AU (Levy & Sonett 
1978). Nevertheless the likely range of field strengths 
is easily estimated. A lower limit is established by the 
1-10 mG magnetic field strengths measured using Zee- 
man splitting of OH 18 cm lines in the cores of molec- 
ular clouds and and in masers in star forming regions. 
This field is likely to be amplified by compression and 
shear during the formation of the central star and its 
surrounding disk. It therefore seems likely that 



B > 10 mG 



(10) 



in protoplanetary disks. An upper limit on B is pro- 
vided by equality of magnetic and thermal pressure in 
the disk midplane, 



87r 



PO Cs 



which yields 



B < B,, 



18 r 



'13/8 
AU 



G. 



(11) 



(12) 



Note that the equipartition field B^q scales as E^/^ if 
the disk surface density is increased over the minimum 
solar nebula value given byeq. Q. The upper and lower 
limits are disparate at 1 AU but are equal at 100 AU. 

It is worth comparing these limits to the field 
strength that is required if magnetic fields are re- 
sponsible for the loss of angular momentum implied 



by the observationally inferred disk accretion rates 
M - 10^"^ Mq yl"^ (Calvet, Hartmann & Strom 2000). 
To do this, we take as a starting point the MHD mo- 
mentum equation 



9v , 

p-^ + p('^- V)v 



VP + pV$ = - 



(B ■ V)B 

"Stt ' 
(13) 



where P is the gas pressure and $ is the point-mass 
gravitational potential of the central star. Assuming 
axisymmetry, and adopting cylindrical coordinates, the 
toroidal component of the momentum equation reduces 
to p [{v ■ \7)v]^ = {B ■ \/B)^ /47r. It is clear that the 
options for angular momentum transport under these 
assumptions are limited! Assuming that Vz is small 
within the disk we obtain 



pvr d{rv^) _ Br d{rB^) B^ dB^ 



dr 



Airr dr 



At: dz 



(14) 



which equates the rate of angular momentum loss by 
the gas to the magnetic torque. The magnetic field 
transports angular momentum radially and vertically 
via the d/dr and d/dz terms on the right hand side 
respectively. The LHS is related to the accretion rate 
at r, 



M = -27rr 



pVr dz . 



(15) 



by integrating between the lower and upper disk sur- 
faces at z — ±s. Before integrating eq ( [T4|) in this 
manner, note that ~ vk4' where vk — \/ GM/r and 
use V • B = to obtain 



pVrVK_ 

2r 



Then integrating vertically between the disc surfaces 
yields 



^'-^§-,i-'h<-BrB, >)-{BM, 
where I have defined 



K—BrBj, >- 



1 

2h 



BrBfjj dz . 



(17) 



(18) 



and BzB^ — +{BzB^)s or —{BzB^)s at the upper or 
lower disc surfaces respectively. If Br and B^ are the 
same order of magnitude, and the magnetic fiux exiting 
the disk surfaces is anchored in an external envelope or 
outfiow, then B^ at the disk surface is non-negligible 



and the first term in (18) is smaller than the second 
by a factor ^ h/r. In this case the torque per unit 
area applied to the disk by the net magnetic tension 
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across it equals the rate of angular momentum loss by 
the material in the disk, and 



<j,\s 



Mvk 
2r2 



1/2 



Q.2M%\-JJ^G, (19) 



where M_7 is M in units of 10""^ Mq yr^-'^. 

This estimate is a robust lower limit on the rms field 
in the disk, as it assumes efficient vertical transport of 
angular momentum by an ordered magnetic field. If ra- 
dial transport dominates instead, then <—BrB(h >^/^ 



must be larger than the estimate ( 19 1 by a factor 
~ \prjh. If the field is disordered because of turbu- 
lence, then these represent a space and time-averaged 
rms magnetic field (see, e.g. Balbus 2003), and the lo- 
cal field may be somewhat larger. For example, simula- 
tions of the magnetorotational instability indicate that 
B„ns 2 <-BrB^ >i/2 (Sano et al. 2004). In any 
case, one concludes that Gauss-strength fields at 1 AU 
are required to give the inferred accretion rates. 



4 Magnetic diffusion 

Magnetic diffusion is an essential ingredient in any the- 
ory of magnetised astrophysical disks. For example, 
field diffusion allows matter to be accreted while leav- 
ing the magnetic field behind. In disk-driven wind mod- 
els magnetic diffusion allows some of the disk material 
to be loaded onto field lines and fiung outwards from 
the disk surfaces (Wardle & Konigl 1993) while the re- 
mainder is accreted. Diffusion is also associated with 
energy dissipation and sets the inner scale of magnetic 
structures. 

The magnitude of the diffusivity is determined by 
microphysical processes rather than the wishes of a 
theorist (or an observer for that matter!). In fully- 
ionised disks the molecular diffusivity is orders of mag- 
nitude too small and an effective turbulent viscosity 
or anomalous diffusivity of plasma physics origin has 
to be invoked as the origin of the breakdown of ideal 
MHD. In the weakly-ionised environment of protoplan- 
etary disks collisions between charged particles and the 
dominant neutrals are sufficient to provide the neces- 
sary diffusivity, and may be too efficient, suppressing 
magnetic activity, leading to the formation of "dead" 
zones in the disk (Gammie 1996; Wardle 1997). If sim- 
ple neutral-charged particle collisions are the origin of 
diffusivity then quantitative calculations can provide a 
sound physical basis for calculations of magnetic field 
evolution in protoplanetary disks. 

The finite conductivity of a weakly- ionised gas can be 
determined by calculating the drift of charged particles 



in response to an applied electric field E' in the neutral 
frame, and summing over the charged species to obtain 
the current (Cowling 1976). The drift velocity Vj of 
each charged species j (mass nij, charge Zjc) relative 
to the neutrals (mean mass m, density p) is determined 
by balancing the force applied by the electric field, the 
magnetic force, and the drag associated with neutral 
collisions: 



Z^eE' + ZjC ^ X B - mj-yjpVj = 



(20) 



where —< av >_,- /{rrij + m) and < av >j is the rate 
coefficient for coUisional momentum transfer between 
species j and the neutrals, so jjp is the collision fre- 
quency with neutrals. Eq (20) implicitly assumes that 



the electromagnetic field and the neutral fluid evolve on 
a time scale that is long compared to the inertial time 
scales of charged particles, which is generally a good 
approximation (see Wardle & Ng 1999). The direction 
of the drift of a particle of mass rrij, charge ZjC is deter- 
mined by the relative magnitude of the magnetic and 



drag terms in eq ( 20 1 , which is characterised by the Hall 



parameter for species j 

^ _ \Zj\eB 1 



(21) 



i.e. the ratio of the gyrofrequency and neutral collision 
frequencjQ Note that apart from molecular factors, the 
Hall parameter depends only on iJ/nn. If (3j ^ 1 the 
Lorentz force dominates the neutral drag 



Z,eE' 



-Z-jC ■ 



X B 



(22) 



and the charged particles are tied to the magnetic field 
lines. In the other limit f3j <^ 1 the drag dominates 



Z,eE' 



■JjfUjpVj 



(23) 



and neutral collisions completely decouple the particles 
from the magnetic field. 



Inverting eq ( 20 1 for Vj , and forming the current 
density J = J2j ''^j^^j'^j yields 



J = aoE'u + ghB X E'^+ apE'^ 



(24) 



where E'^^ and E'_^ are the components of E' parallel 
and perpendicular to B and the ohmic. Hall, and Ped- 
ersen conductivities are (e.g. Cowling 1976; Wardle & 
Ng 1999): 



(25) 



^Note that Wardle & Ng (1999) include the sign of Zj in the 
definition of /Bj 
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and 



(Tp 



B 



ec 
~B 



St 



B 



-/5 



(26) 



(27) 



respectively, where I have used ^ rij Zj = in deriving 



the second form of (Jh in eq ( 26 1 



Inverting eq ( 26 1 for E' yields an induction equation 
of the form 



dB 



= W X {v X B)-W X [770V X B 
+ 77h(V X B) X B + 77a(V X B)j_ 



where 



Vo 



47r (To 



and 



47rCT_L crj_ ' 



ap 



Vo 



(28) 



(29) 



(30) 



(31) 



4tT(T± (Tj_ 

are the ohmic, Hall and ambipolar difFusivities respec- 
tively, where a± — \J a\ + a\, . 

If the only charged particles are ions and electrons, 
then 77H = /3e'7o and ryA = PiPeVo- The Hall parameters 
for ions and electrons are 



P, « 4.6 X 10' 



and 



/3e « 3.5 



n-15 



T 



(32) 



i3G 
ni5 VlOOK 



-1/2 



respectively, where is the magnetic field in Gauss 
and ni5 = nn/lO^^cm"^. As (3elPi ~ 1000 there are 
three distinct diffusion regimes: 

A ^ /?e ^ 1 ohmic (resistive) 
A < 1< /3e Hall (34) 
1 (3i Pe ambipolar 

which correspond to different regimes of B/nn, as il- 
lustrated in Figure [T] In the ohmic regime, the ion 
and electron drifts are unaffected by the magnetic field. 



en 
o 




log n 

Fig. 1. — The magnetic diffusion regimes of a weakly- 
ionised, three-component plasma are determined by the 
ion and electron Hall parameters /3i and /3e, which are 
proportional to B/n, with Pe/Pi typically ^ 1000 (see 
text). 

in the Hall regime the electrons are tied to the field 
but the ions are not, while in the ambipolar diffusion 
regime both species are tied to the magnetic field and 
drift together through the neutrals. The presence of 
grains complicates this picture (Wardle & Ng 1999), 
but generically ohmic or ambipolar diffusion dominates 
if the majority of charged particles are tied to the neu- 
trals or the magnetic field respectively, otherwise Hall 
diffusion is important. 

Which diffusion regimes are relevant? Figure |2] plots 
the loci in the log _B-log nn plane where the Hall param- 
eters are unity for ions, electrons, and singly-charged 
grains with radii of 50 and 2500 A. Also shown are 
the regions occupied by molecular clouds (_B[mG] ~ 
•^/riHicm"-^] ) and the midplane of the solar nebula (eqs 



and 12 ) between 1 and 100 AU. Densities are so 



large at the midplane of the solar nebula that ions are 
decoupled from the magnetic field by neutral collisions. 
This is not so for electrons except for very weak fields, 
and so Hall diffusion is important once grains have set- 
tled or aggregated and are not the dominant charge 
carriers. If grains are important they are very much 

(33) 

decoupled from the magnetic field (i.e. Pg <C 1) and 
ohmic diffusion dominates. However, the diffusivity is 
then so severe that the magnetic field cannot couple 
effectively to the gas at all (Hayashi 1981). 

By contrast, in molecular clouds ions and electrons 
are tied to magnetic fields whereas the largest grains are 
not. The charge residing on large grains is usually small 
except inside shock waves (Wardle & Chapman 2006) 
or at very high densities. Thus in molecular clouds, 
ambipolar diffusion dominates. 

How much diffusion in protoplanetary disks is too 
much of a good thing? This depends on the scale on 
which one would like the magnetic field to couple to 
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log Hh (cm ^) 

Fig. 2. — Loci in the log B-log nn plane where the Hall 
parameters are unity for ions and electrons (red) and 
50 and 2500 A grains (blue). Thick black lines indi- 
cate the region occupied by molecular clouds (-B[mG] ~ 
y^TjflJcrrr^) and the equipartition field at the mid- 
plane of the minimum- mass solar nebula between 1 and 
100 AU. 



the matter in the disk. The weakest criterion for in- 
teresting magnetic effects is to demand that the mag- 
netic field should at least be able to couple to the Ke- 
plerian shear in the disk. This is required both by 
the magnetorotational instability and by disk-driven 
wind models. To estimate this, consider the induction 
equation ( 28 1 . The magnitude of the inductive term 
|V X (v X B)| ~ QB, where fl = vk/t = Cg/h is the 
Keplerian frequency. The diffusive term is dominated 
by vertical gradients, and so |V X {t]\/ X B) \ ^ -qB/h? 
where 77 is the largest of 'qo, \'<1h \ and rjA- The diffusion 
term must be smaller than the advective term to allow 
the field to couple to the Keplerian shear, correspond- 
ing to the requirement 



77 < hcs (35) 
for good coupling. 

5 Diffusivity calculations 

To answer the question of whether the magnetic diffu- 
sivity is too low or high, and which diffusion mechanism 
dominates in different regions of protoplanetary disks, 
we must calculate the ohmic. Hall and ambipolar dif- 
fusivities appearing in eq ( 28 1 . These depend on the 



abundances of the charged species and their Hall pa- 
rameters (see eqs 



29 - 31 ) 



is isothermal in the z-direction. The ionising sources 
are radioactivity (which is largely negligible), cosmic 
rays, and the x-ray flux from the central star. The cos- 
mic ray flux is assumed to induce an ionisation rate 
10~^^s~^H~^ well outside the disk and suffers expo- 
nential attenuation with depth as exp(— E/96gcm~^) 
(Umebayashi & Nakano 1981). For the depth depen- 
dence of the x-ray ionisation we adopt a fit to the results 
of Igea & Glassgold (1999). The simple chemical reac- 
tion scheme of Nishi, Nakano & Umebayashi (1991) and 
Sano et al. (2000) is adopted, which follows the abun- 
dances of the light ions H+, H^, He+, C"*", and repre- 
sentative heavier molecular and metal ions (denoted m+ 
and M"*" respectively), as well as charged grains. Their 
scheme has been extended to include high grain charge 
using the cross sections of Draine & Sutin (1987). This 
is necessary because the temperatures here are high 
compared to molecular clouds - the high thermal veloc- 
ities of electrons can overcome large Coulomb barriers 
allowing grains to pick up a high charge through elec- 
tron sticking. Following Umebayashi & Nakano (1990) 
I consider populations of uniform size grains, as well as 
models in which grains have settled out. 

A sophisticated treatment of the ionisation equilib- 
rium using a full chemical reaction network has been 
performed by Semenov, Wiebe & Henning (2004), al- 
though with a more restricted treatment of grain charg- 
ing. 




To calculate the charged particle abundances at a 
given radius we adopt the surface density and temper- 
ature of the minimum solar nebula, and assume that it 



Fig. 3. — Ionisation rate (black dotted curve) and frac- 
tional abundance (relative to hydrogen nuclei) of ions 
(blue) and electrons (red) as a function of height above 
the midplane, z, in the minimum solar nebula at 1 AU 
from the central star. Grains are assumed to have ag- 
gregated and settled to the midplane. 

Figure [3] shows the ionisation rate and abundances 
at 1 AU as a function of height z from the midplane in 
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the absence of grains (i.e. once they have settled). The 
ionisation rate is dominated by x-rays above z/h — 2, 
and by cosmic rays at greater depths. The dominant 
charged species are electrons and metal ions because 
the recombination rate coefficient for metal ions is two 
orders of magnitude below that for molecular ions. The 
ionisation fraction declines with depth because of the 
decreasing ionisation rate and increasing neutral den- 
sity. 

Figure |4] shows the corresponding ohmic, Hall and 
ambipolar diffusivities as a function of height for field 
strengths of 0.1 and 1 G. The ohmic resistivity is inde- 
pendent of the magnetic field strength, whereas the Hall 
and ambipolar diffusivities scale linearly and quadrati- 
cally with B respectively. The ohmic diffusivity is pro- 
portional to Tig and declines strongly with height. The 
ion and electron Hall parameters, which are inversely 
proportional to p, increase strongly with height. As 
a result, the Hall and ambipolar diffusivities increase 
above z/h ^ 2 where the neutral density starts to 
drop rapidly. Hall diffusion dominates below 2-3 scale 
heights, and ambipolar diffusion dominates at greater 
heights. Note that the magnetic field is able to couple 
to the shear in the disk below 3-4 scale heights. Above 
this, ambipolar diffusion is severe. At first sight this 
is counter-intuitive - the diffusion is most severe above 
the disk where the ionising flux is strongest and the 
fractional ionisation is greatest! The problem is that 
the number density of charged particles is very low at 
the disk surface. 

Figure[5]extends this plot to other values of B by rep- 
resenting the diffusivity as contours of ('7o + ^h+^a)^^^ 
in a log B-z/h plane, with the background colour indi- 
cating the dominant diffusion mechanism. The curves 
in Figure [4] correspond to horizontal cuts in this plane 
at log i?(G) — and —1. It is clear that ohmic diffu- 
sion is important only for weak fields near the midplane 
(i.e. when f3e ^ 1), whereas ambipolar diffusion is dom- 
inant at the disk surface and for stronger fields (when 
/3i ^ 1). Hall diffusion dominates over the large range 
of intermediate conditions. 

If we add a standard interstellar population of 0.1/im 
radius grains, the ionisation equilibrium is very different 
(Fig. |6| because the grains acquire a charge via stick- 
ing of electrons and ions from the gas phase. Above 
z/h K A the grain charge is determined by the compet- 
itive rates of sticking of ions and electrons. The mean 
grain charge, « — 7e, is determined by the Coulomb 
repulsion of electrons offsetting their greater thermal 
velocity compared to ions. Most recombinations still 
occur in the gas phase. For z/h « 3-3.5, the abun- 
dances of ions and electrons have declined to the point 
that the majority of electrons stick to grain surfaces be- 




z/h 

Fig. 4. — Magnetic diffusivities (ohmic - black curve; 
Hall - red curves; ambipolar - blue curves) as a func- 
tion of height above the midplane derived from the 
charged particle abundances plotted in Fig. [3] for mag- 
netic field strengths of 0.1 and IG (solid and dashed 
curves respectively). A single curve for ohmic diffu- 
sivity is plotted as this does not depend on B. The 
horizontal dashed line indicates the value of diffusiv- 
ity at which the magnetic field is unable to effectively 
couple to the Keplerian shear in the disk (see ^4|. 

fore they can recombine in the gas phase, and most neu- 
tralisations occur when ions stick to negatively charged 
grains. Closer to the midplane, the ionisation fraction 
is so low that most grains are neutral, and ions typi- 
cally stick to a neutral grain before encountering a neg- 
atively charged grain. Recombinations mostly occur in 
this region through the collision of positive and nega- 
tive grains. These regimes were found by Umebayashi & 
Nakano (1990), and as expected the curves are similar 
to theirs (except for the large grain charge for z/h> 3) 
when the abundances are plotted as a function of C/"-h 
rather than z/h. 

Grains severely decrease the conductivity by reduc- 
ing the abundance of charged particles and because of 
their much reduced mobility relative to ions and es- 
pecially electrons. Figure [7] shows the diffusivity as a 
function of height for field strengths of 0.1 and 1 G. 
Close to the midplane the diffusivity is increased by 
7-8 orders of magnitude compared to the case when 
grains are absent. It is only for the weaker field and 
z/h > 3.5 that the diffusivity is low enough to allow 
the magnetic field to couple to the disk. The diffusiv- 
ity as a function of z/h and field strength is plotted in 
Figure [8] The vertical contours on the left hand half of 
the plane (i.e. for z/h < 3 refiect the strong decline of 
the conductivity towards the midplane as the density 
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Z/h 



Fig. 5. — Logarithmically-spaced contours of (77Q + 
+ '7a)^^^ ^ ^ function of height and magnetic field 
strength at 1 AU in the minimum mass solar nebula. 
The charged particle abundances are plotted in Figure 
|3] The contour levels increase by factors of 10 from 
lO""*^ hcs at the bottom of the plot near = 3 up to 
10^ hCs at the very top right-hand corner of the plot. 
The solid contour is the critical value hcs - in the re- 
gion above this contour the magnetic field cannot cou- 
ple effectively to the disk. The background shading 
indicates whether the dominant diffusion mechanism is 
ohmic (purple), Hall (blue) or ambipolar (green). 
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Fig. 6. — As for Fig. |3] but including a standard inter- 
stellar population of 0.1 /im radius grains. Green curves 
indicate the abundances of different grain charge states 
Ze, labelled by Z. 

(and therefore neutral drag) increases and the num- 
ber of charge carriers decreases. For weak to moderate 
fields the dominant charge carriers - i.e the grains ~ are 




Z/h 

Fig. 7. — As for Fig. |4] but for the charged particle 
abundances plotted in Fig. |6] Note that the vertical 
scale extends to much larger values than that in Fig. |4] 
because the dominance of charged grains for z/h < 3.5 
increases the diffusivity by many orders of magnitude. 
The kink in the dashed red curve at z/h w 2.2 occurs 
because for S = 1 G the Hall diffusivity becomes nega- 
tive above this point (the absolute value is plotted). 

decoupled from the magnetic field by neutral collisions 
and ohmic diffusion dominates. The grains become tied 
to the magnetic field if it is strong enough and ambipo- 
lar diffusion dominates. Hall diffusion, which would 
normally be prominent at intermediate field strengths, 
is suppressed by the almost equal numbers of positive 
and negative grains. This is because the Hall effect 
relies on different degrees of field line-tying for posi- 
tive and negatively charged species. Above z/h k, 3.5 
where grains become unimportant, the diffusivity be- 
haves as in Figure [5] The region where the magnetic 
field can couple to the differential rotation of the disk 
is restricted to B < 0.5 G and z/h « 3.5-5. 

These effects can be mitigated by the aggregation 
of grains, which reduce the surface area available to 
capture charged particles from the gas phase. For ex- 
ample, I show in Figure [9] the diffusivities obtained if 
grains have aggregated to 3/Ltm in radius. This is suf- 
ficient to allow the coupled region to extend down to 
z/h « 1.5, which is now the height below which the 
grains dominate the charged species. Grains rapidly 
grow to this size in protoplanetary disks, so this sug- 
gests that Gauss-strength magnetic fields can couple to 
10% of the disk at 1 AU. 

At 5AU, the region of giant planet formation, the 
x-ray flux is reduced by a factor of about 25, while 
the column density, « GOgcm^'^, is no longer able 
to significantly attenuate the cosmic ray flux imping- 
ing on the midplane. The ionisation rate drops from 
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1 AU 0.1 Mm grains 




5AU 



no grains 



Fig. 8. — As for Fig. |5] but for the charged particle 
abundances plotted in Fig. [6] (i.e. including 0.1 /im ra- 
dius grains.) 




Fig. 9. — As for Fig. [5] but assuming that grains have 
radius 3 /im. 



the midplane. The ionisation fraction is similar to that 
at 1 AU, but the reduced gas density allows greater mo- 
bility of ions and electrons. Equipartition fields 1 G) 
can easily couple to the midplane in the absence of 



grains (see Fig. 10). Note that Hall diffusion tends to 
dominate. 

Grains, if present, increase the diffusion and create 
a dead zone, though not as extensive as that at 1 AU. 
Increasing the grain size to just 1 /im is sufficient to 
move the transition from gas-phase to grain-dominated 
recombinations down to the midplane. Then magnetic 
fields with strengths below a few tens of milliGauss can 
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Fig. 10. — As for Fig. |5] but at 5 AU in the minimum 
solar nebula in the absence of grains. 



couple to the bulk of the disk (Fig. 11) and Hall diffu- 
sion again is dominant. 
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Fig. 11. — As for Fig. [5] but at 5 AU in the minimum 
solar nebula and assuming that grains have radius 1 fim. 



6 Discussion 

These resistivity calculations assume a minimum-mass 
solar nebula subject to cosmic-ray and x-ray ionisation, 
with grains (if present) assumed to be of uniform size. 
Here I discuss the effect of relaxing these assumptions. 

First, it is not at all clear that protoplanetary disks 
are subject to ionisation by interstellar cosmic rays be- 
cause super- Alfvenic winds may advect them away from 
the disk (as noted by Gammie 1996 and Glassgold et 
al. 1997), or they may be excluded by unfavourable 
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magnetic geometry such as hour-glass shaped magnetic 
field lines vertically threading the disk or a wound-up 
field lying in the disk plane. This may be compen- 
sated by cosmic rays produced by magnetic activity 
near the central star or in a disk corona, but it is also 
unclear whether these particles are able to impinge on 
the disk. If cosmic rays are absent then the only signif- 
icant ionisation source are x-rays from the central star. 
Naively one might think that the active layer is then 
restricted to the nominal penetration depth of x-rays, 
i.e. ~ lAmji/aT ^ 3.5 gcm^^ . However, the x-ray 
induced active layer is much deeper than this because 
the unattenuated x-ray ionisation rate is so high - x- 
ray ionisation dominates for z/h > 2. I find that in 
the absence of dust grains, the active column density is 
Sactivo ~ 150gcm~^. When dust grains are present the 
active layer lies within the region where x rays domi- 
nate cosmic rays in any case, so this case is unaffected 
by eliminating the cosmic rays. 

Second, the calculations presented here adopt the 
surface density for the minimum-mass solar nebula. 
While this is broadly consistent with observed proto- 
planetary disks (Kitamura et al. 2002), it is possible to 
have surface densities 30-50 times higher without suf- 
fering gravitational instability. The approximate effect 
of adding material to the disk is to multiply up the 
density everywhere by the same factor. The dominant 
factor in the ionisation equilibrium is the shielding of 
the ionising fluxes of cosmic and x-rays by the overlying 
column of disk material. To a first approximation the 
active column densities discussed above are unchanged. 

Third, I have implicitly used single size grain models 
to crudely mimic the effect of aggregation. Some care is 
needed in interpreting these results because of the ten- 
dency for the smallest grains to carry a significant frac- 
tion of the grain charge. For example, although the typ- 
ical grain size increases to micron scales on times that 
are short compared to the disk evolutionary time scale, 
the residual small grains implied by any stochastic ag- 
gregation process may still be important (see, e.g. No- 
mura & Nakagawa 2006). The primary role of grains is 
to mop up electrons from the gas phase and render them 
immobile, increasing the magnetic diffusivity by many 
orders of magnitude. This effect is significant below the 
transition layer where the charge density of electrons 
and grains are comparable (e.g. at z/h « 3.5 in Fig. [6]). 
Above this transition, the charge residing on grains is 
unimportant in determining the abundances of ions and 
electrons and their occasional sticking onto grains yields 
a mean charge Zg{a) e « —AkTa/e on grains of radius a 
(Spitzer 1941; Draine & Sutin 1987). The total charge 
per unit volume carried by grains with a size distribu- 
tion Xg{a) (where the number density of grains with 



radii between a and a -I- da is nH-Tg(a) da) is therefore 
Jxg{a)Zg{a) da. The transition from low to high diffu- 
sivity occurs when the charge per hydrogen nucleus on 
grains becomes comparable to that in electrons, so the 
height of the magnetically active surface layer (zmagi 
say) is a weak, monotonically increasing function of 
JaXg{a) da. The results presented for the 3 fim and 
0.1 /im grain models show that Zmag/^ ^ 2.2 and 3.5 
for Jaxg{a)da = 8.9 x lO^^^cm and 9.9 x 10"^^ cm 
respectively. By way of comparison, 1% of the total 
grain mass residing in O.Ol/xm grains would carry ap- 
proximately the same net charge as the 0.1/xm grain 
population so a trace residual population of small grains 
may poison the magnetically coupled layer. 

Finally, although the fractional ionisation increases 
with height above the disk the diffusivity also increases. 
This is because the number density of charged par- 
ticles plummets because the density is dropping as 
exp(— 2:^/2/1^). This potentially affects the ability to 
drive a wind from the disk surface. However it should 
be noted that the density profile at large heights be- 
comes increasingly unreliable due to the effects of tem- 
perature stratification and/or the presence of outfiows. 

7 Summary 

If dust grains have settled to the midplane, or have oth- 
erwise been removed from the gas, the ionisation level 
in the disk is sufficient to couple the magnetic field can 
to the midplane even within 1 AU of the central star. In 
this case the column density of the magnetically active 
region of the disk is Sactivc ~ 1700 gcm^^. Hall diffu- 
sion dominates (Wardle & Ng 1999; Sano & Stone 2002) 
because the ions are decoupled form the magnetic field 
and electrons are not. Notably, the Hall effect imparts 
an implicit handedness to the fluid dynamics, which be- 
comes sensitive to a global reversal of the magnetic fleld 
direction. 

Dust grains, if present, substantially increase mag- 
netic diffusion by soaking up electrons and ions from 
the gas phase. Grains have huge cross sections for neu- 
tral collisions and this combined with their large inertia 
effectively decouples them from the magnetic fleld and 
the fluid becomes resistive. Magnetic diffusion is so 
severe that the fleld is unable to couple to the shear 
in the disk except close to the disk surface. For a 
standard interstellar population of O.l/im-radius grains 
the active surface layers have a combined total column 
5^active ~ 2gcm^^. As grains aggregate, their effect on 
the ionisation equilibrium is reduced, and by the time 
grains have aggregated to 3 ^m, Sactivc ~ 80 g cm~^. At 
5 AU the lower column density and reduced gas density 
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and shear means that the coupUng is more easily main- 
tained. In this case Sactive ~ 5]totai once grains have 
aggregated to 1 /xm in size. 

If a substantial grain population is present, ioni- 
sation in the magnetically active layers is dominated 
by x-rays so that there is little change if cosmic rays 
are unable to impinge on the disk. In the absence of 
grains and cosmic rays, stellar x-rays maintain mag- 
netic coupling to 10% of the disk material at 1 AU (i.e. 
^active ~ 150gcm~^) and almost all of it at 5 AU once 
grains have aggregated or been removed by settling, 
which is expected to happen on a time scale of a few 
thousand years (e.g. Nomura & Nakagawa 2006). 
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